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Abstract 

Statistical mechanics of a disordered system of cars on a single-lane road is de- 
veloped. Behaviour of cars is defined by conditional probability of car velocity 
depending on the distance and velocity of the car ahead. A system consisting of 
different cars is modelled by a system of two types of cars differing in maximal 
velocity or efficiency of brakes. Starting from conditional probabilities and using 
principle of maximum entropy, probability densities of car velocities and headways 
are calculated. It is shown that the first-order phase transition between free flow 
and congested traffic may be driven by number of fast cars in a system of slow cars, 
and, as a rule, admixture of cars of superior qualities does not increase but decreases 
the total flow. In the system of cars with poor brakes platoons of cars of the same 
velocity are formed. They are dissolved by a small addition of cars with good brakes. 
Application of principle of maximum entropy was justified by comparing the results 
with steady state properties of an equivalent kinetic model. 
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1 Introduction 



The behaviour of cars on a single-lane road without possibility of overtaking 
is given by reaction of each driver on the velocity and distance (headway) of 
the car ahead. The drivers cannot influence the headway directly, they control 
only the velocity which is limited by the maximum construction velocity of 
the cars and necessity to avoid an accident with the car ahead if it is slowing 
down or stops. The driver does not react on the traffic situation always in 
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the same way, so that the velocity of his car would be distributed around 
some optimum value with a probability that can be found experimentally or 
modelled, starting from simple considerations. At the beginning the model is 
formulated as continuous, with arbitrary coordinates and velocities from given 
intervals, nevertheless, the numerical calculations are performed on a one- 
dimensional lattice where the car acquires only discrete velocities. Then the 
model resembles discrete particle hopping models [1,2,3], but the configuration 
space of particles is enlarged from a coordinate to the coordinate and velocity 
space, and instead of hopping probabilities, particle velocities are introduced. 
The model is apphed to a steady state traffic, i.e., to a system of cars of 
constant global density observed for a long period of time. In the real traffic 
the cars are not identical. This fact and its infiuence on the distribution of 
velocities and headways of the cars is exemplified by investigation of a random 
mixture of two kinds of cars. 

Our model of traffic and the method of its solution is described in Section 2. In 
Section 3, the results in form of 3D plots of probability density of car velocity 
and headway averaged over a small groups of cars are presented. 



2 Model of traffic flow 



The behaviour of the i-th car is further described by a conditional probability 
density pi{vi\vi-i,Xi-i^i) of its velocity Vi depending on the distance (head- 
way) and velocity of the preceding car [4]. The probability density pi can be 
found experimentally observing a couple of cars for a long enough period. The 
knowledge of the probability gives us full information of driver's response to 
traffic situations, and no assumptions on reaction times, breaking abilities, 
engine power, etc., are necessary. The driver can directly influence only the 
velocity of the car, not its relative position with respect to its neighbours. 
There is no a priori information about the distances between the cars, and 
they will be obtained as a result of calculations. This fact is taken into ac- 
count by maximization of Shannon's information entropy which is identical 
with Gibbs entropy in statistical mechanics. The probability distribution of 
distances is obtained from it together with the requirement that all the cars 
have the same mean velocity as they cannot overtake each other. 

The states of the system are characterized by the values of velocities and 
coordinates of the vehicles, and our task is to find the probability of each of 
such states, 2^1,2, '?^2! • • • ,xn-i,n,Vn)- (The cars are numbered from the 

left to the right and are moving to the left.) As each car interacts only with 
the preceding car, fixing the velocity of one car divides the system into two 
independent parts, i.e.. 
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V{vi, xi,2, V2,...,Vi^ const, . . . , Xn-1,N, Vn) = 

= 7^1(^1, Xi^2, V2,..., Xi-i^i)V2{Xi,i+l, Vi+i, Xn-1,N, Vn) 



where V2 depends on the value of Vi and is, in fact, a conditional probability. 

Applying the above consideration to each vehicle, we see that the probability 
that the velocities of cars in a group are Vi and the distances between them 
Xi-i^i is a product of conditional probabilities 

N 

V{vi, Xi,2, V2,..., Xn-1,N, Vn) = Pl{vi) n Pi{^i-l,i^ Vi\Vi-i), (l) 

i=2 



where 



PiiVi_i, Xi_i^i, Vi) = Pi{Vi\Vi_i, Xi_i^i)Pi{Vi_i, Xi_i^i), 

Pi{vi-i) = j dxi^i^iPi{vi-i,Xi-i,i) = 

Pi-l{Vi_i) = j <lXi_2,i-lPi-l{Xi^2,i-l,Vi^i). 



(2) 



Substituting for Pi{xi-i^i,Vi\vi-i) from (2) into (1) we obtain 



V{Vi, Xi^2, V2,---, Xn-1,N, Vn) 



N 



Pl{vi) '[[pi{Vi\Vi-i,Xi-i,i) 



Pi{Vi^l,Xi^l^i) 



1=2 



Pi{Vi 



(3) 



The probability V has a relatively simple form as we assume that the driver 
of i-th car reacts only on the velocity of the i — 1-th car and its distance 
and not on the velocities and distances of the other cars. The conditional 
probability density Pi{vi\vi-i,Xi^i^i) contains all the known information of the 
system, but it cannot alone determine the probability of states of a group 
of cars. To find it, having no further information available, the principle of 
maximum entropy is used. To calculate the probability V, we have to maximize 
entropy log V under conditions that the mean length of the system is 

L and the mean velocity of each car is v, which, in equilibrium, has to be the 
same for each car. Using the method of Lagrange multipliers, we have to find 
the extreme of the expression 
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J l[dxi,i+idviV\ogV + 

i 
i 

+ Y^^i j ^^i'^i j / dXi_i^iPi{Vi\Vi_i,Xi_i^i)Pi{Vi_i,Xi_i^i) + (4) 

i 

where n, i/i, Aj are Lagrange multipliers. 

Substituting (3) into (4) and maximizing it with respect to Pi{vi-i,Xi_i^i), we 
get 

log _^ i_ixi_i^i + PiVi{vi-i, Xi-i^i) + Si{vi-i, Xi^i^i) + Ai = (5) 



where Vi(vi^i,Xi^i^i) = J dviViPi{vi\vi^i,Xi^i^i) and Si{vi-i,Xi^i^i) = 
J dviPi{vi\vi-i,Xi-i^i) ■ logpi{vi\vi-i, Xi-i^i) are the mean velocity and entropy 
of a car with fixed boundary conditions. / dxi-i^i Pi{vi-i, Xi-i^i) = Pi-i{vi-i) 
is also required. Then, using (2), (3) and (5), the probabilities appearing in (1) 
can be calculated from the following formula 

Pi{xi-i,i,Vi\vi_i) ^ Pi{vi\vi^i,Xi_i^i)Pi{xi^i^i\vi^i) (6) 
Pi{xi-i,i\vi^i) = exp(-//Xi_i,j - iyiVi{vi^i,Xi^i^i) - Sj(^;i_i, 

where A{vi-i) is the normalization constant 



j dxj_i,iexp(-/xxi_i,, - UiVi{vi^i,Xi^i^i) - Si{vi^i,Xi^i^i)) 



Prescribing the mean velocity and density of cars, the Langrange multipliers Ui 
and II, respectively, can be calculated from these requirements. In a system of 
identical cars all Ui would be the same and equal to global u. As the velocity is 
a nonconserved quantity, z/ = in full equilibrium, where the velocities are not 
given, but obtained from the requirement of maximum entropy. Multipliers i^i 
are non-zero for cars with an obstacle in front of them, causing reduction their 
velocity, or if the cars are not identical, and the slower cars are hindering the 
faster ones. 

The statistical properties of a car are given by the probability density that 
the car velocity is Vi and its headway Xi^i^i 

Pi{xi-i,i, Vi) = j dvi^i Pi{xi^i,i, Vi\vi_i)Pi^i{vi_i) (7) 
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where Pi{xi-i^i,Vi\vi-i) is calculated from (6). 

The probability density Pi{xi_i^i,Vi) reflects the behaviour of a single car. 
Correlations between velocities and positions of different cars are given by 
two- or many-body correlation functions. We shall not calculate them here. 
Instead of it, the probability density P{V,X) that the average car velocity 
in a small group of m cars is V = (1/m) X]™ i f j and average headway X = 
(1/m) J2iLi is calculated [6]. This quantity yields the information about 
the distribution of single car properties as well as correlations between the cars. 
The contribution of correlations to the distribution is given by the difference 
between the probability density P(V,X) and the probability that a group of 
uncorellated cars with probabilities Pi{xi-i^i,Vi) have the same values of the 
averaged velocity and headway 



N 

D{V, X) = jU dxi_i,dvi S{V - J2 v^/m)6{X - ^ x,_,,/m) 

i=l i i 

.AT N 

/ dvo'[[Pi{xi_i^i,Vi\vi_i)P{vo) - '[[Pi{xi-i^i,Vi) 

i=l i=l 



(8) 



The conditional probability of velocity of the i-th car for given headway and 
velocity of the preceding car Pi{vi\vi_i, .Xj_i j) is not acquired from experiment 
in this paper, but calculated from simple considerations [4] , close to the ideas 
of car- following model [5] . It is assumed to be peaked around an optimal ve- 
locity Vi^opt, which is further chosen as 90% of maximum safe velocity Vi^rn- The 
maximum safe velocity is determined from the requirement that two neigh- 
bouring cars, which start to decelerate at the same time with their deceleration 
rates aj_i,aj, would stop without crash. Moreover, Vi^^a must not be greater 
than the maximum possible velocity of the car Vi,max, i-e., the optimal velocity 
of car i is 



Vi,opt{vi-i,Xi-i,i) = 0.9vi^,^{vi^i,Xi^i^i), (9) 



-a-iTi \j{aiT^'^ 2aiXi_x,i + ^'i^-iOi/ai-i if Vi,™ < ^"1,11 

^i,max i-f ^i,m ^ ^i,E 



where a;j_i^j and are the distance (headway) and velocity of the car ahead, 
respectively, Oj denotes the deceleration rate of car j in time of braking, and 
Tj is the reaction time of the j'-th driver. 

The cars in (9) are not identical. They differ in their maximum velocities and 
deceleration rates. In the formula, it is assumed that the driver of a car can 
distinguish the type of the preceding car and knows its maximum deceleration 
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rate. In the opposite case the driver i should expect that Oj-i = max aj for all 
J- 

The equations of the above-described model are further solved numerically. 
For numerical computations, it is more convenient to treat the model on a 
discrete lattice with velocities acquiring only integer values. This approach 
was also used in kinetic models of traffic flow [1]. Now, our model may be 
considered as a generalized version of the ASEP model [7] . It was shown that 
our equihbrium model describes well the steady states of the ASEP model. 
The space dependence of density and pair correlation functions is very close 
to the exact solutions [4] . 

In numerical calculations, the cars are presented by dimensionless points mov- 
ing on a discrete one-dimensional lattice and are characterized by 2 quantities: 
discrete integer velocity Vi in the interval (0, Vj^max) and a discrete coordinate 
(site number) Xi = 1,2,... (further, only their differences will be con- 

sidered). fj,max is the maximum velocity given by the construction of the car 
(the nearest integer value to Vopt from (9) is taken for the optimal velocity). 

As the road is assumed to be homogenous, in the calculation no absolute 
coordinates are used, only the distances between the cars. Nevertheless, the 
velocities are absolute, measured with respect to the road. 

If a homogeneous system is in a steady state for a long time, the small groups of 
m cars may be chosen at arbitrary places and arbitrary moments. In practice, 
we would observe small parts of the road (e.g. by TV cameras) for the time 
the density of the whole system remains constant. 

The way of driving of a driver is characterized by distribution of probabili- 
ties of car velocities around the optimal velocity. Here we use an extremely 
simple distribution, in which the probability of the optimal velocity is po, the 
probabilities of the velocities fopt =t 1 are pi, while the probability of the car 
to have any other permitted velocity is p2. The values of the probabilities for 
velocities higher than the maximum safe velocity are equal to 0. The sum of 
all probabilities for each car is equal to 1. The parameters po,pi, and p2 are 
the same for each car, and the distribution depends on the headway only by 
means of the value of optimal velocity. For moderate densities the continu- 
ous model is described well by our discrete model defined by the probabilities 
p. Nevertheless, for very high densities, the discreteness of the lattice is not 
negligible, and another condition Vi < Vi-i + Xi-i^i must be imposed. 

Our approach enables us to treat a disordered system with each car having 
different properties. To exemplify the effect of disorder, we shall investigate a 
system consisting of two types of cars for two simple cases: cars are differing 
by i) maximum velocity, ii) efficiency of brakes. 
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In an equilibrium state of a large group of cars, the cars far from the foremost 
edge of the group have the same mean velocity. It is calculated by an iteration 
starting from the first car of the group driving, e.g., by constant, optimum 
velocity ^1=^1 = fopt, i.e., -Pi(fi) = Svi,vopf Then, the velocity distribution 
of the second car is P2{v2) = T,vi,xi,2 P2{xi,2:V2\vi)Pi{vi) where P2{xi^2:V2\vi) 
is given by (6) with z/2 = if ^2 = E^a '^2-P2('y2) < '^i- If ^2 > Vi, P{v2) is 
recalculated with 1/2 > slowing down car 2 to the velocity of the preceding 
car (1), i.e., V2 = Vi. This ^(^2) is used in the next step. The same procedure 
is applied consequently to all other couples of cars i — l,i in the group. After 
large enough number of steps the mean car velocity converges to a constant 
value Vi = V. The density of cars is controlled by the Lagrange multiplier 
(potential) ^ which is the only free parameter of the theory (if the conditional 
probabilities Pi are known from the above considerations or an experiment 
with two cars), i.e., also v is a. function of /i. The Lagrange multipliers are 
determined from the requirement Vi — v. For a system consisting of identical 
cars, all z/j = if the velocity of the first car is large enough. In a mixed system 
PiS of the fast cars are positive to slow them down to the mean velocity of the 
cars in front of them. If there are no obstacles slowing down the velocity of the 
group, at least some multipliers Ui of the slow cars are equal to zero. In the 
further considerations, the mean car velocity is not confined, and its natural 
value assuming free road without any velocity limitations, except those due 
to construction of cars, is calculated. 

After reaching the steady state, probability distribution P{V,X) of average 
velocity F of a car and average headway X in a small group is calculated. It 
corresponds to an experimental observation of a constant density traffic with a 
TV camera displaying only a small part of the road. Moreover, for disordered 
systems with two kinds of cars, the probability distribution is averaged over a 
large number of randomly selected small groups. 

To justify the use of the principle of maximum entropy in the equilibrium state 
of our system of cars, the results were compared with the results obtained 
from numerical simulations of steady state of an equivalent stochastic kinetic 
model. In this model particles occupy sites of a discrete ID lattice and perform 
jumps in one direction of length Axi{t) at discrete times t with conditional 
probability p,;(A:r,;(/:)| A:r,;_i(f:), Xj — 1)), where Xi^i^i{t — 1) = — 
2) + Axi(t — 1) — Axi^iit — 1) foT i = 2, . . . , N. The function pi is identical to 
the conditional probability of the velocity of the i-th car, pi, defined above and 
used further in calculation of equilibrium properties. The jumping probability 
of the first car in the group, Pi{Axi{t)) is calculated selfconsistently as a 
probability of the jumping length of the cars of the same kind in the group 
of cars at the time t — 1. The probability densities of the headway X 
and car velocity V were averaged over the last five cars of a large group 
and large enough numbers of time steps after reaching the steady state. The 
probability densities obtained from the steady state of the kinetic model and 
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from the equilibrium distribution assuming the principle of maximum entropy 
were compared for both ordered and disordered systems and found to be very 
close to each other. 



3 Results 

The system is described by a large number of parameters. To study only 
few typical situations some of them are further fixed: pi/po = 0.3, r = 0, 
and the number of cars in the group, m — 5. The others: spread of the car 
velocities P2/P0, its maximum velocity v^ax, deceleration rate a, potential /i, 
and concentration ratio of the two types of cars will acquire a small number 
of values. Typical traffic situations, shown in the following plots, appear at 
Wmax = 12, 20, = 0.07, 0.1, 0.12, a = 0.5, 4, and ps/po = 0.33, 0.55. 

The calculations performed for systems consisting of identical cars show [4] 
that they might be found in a free flow regime or a congested phase. For 
very high densities, platoons of cars of the same velocity are formed. The free 
flow and congested phase may coexist, and a first order phase transition takes 
place. It is reflected in the shape of the probabihty density P{V, X) with two 
peaks. The relative height of them is controlled by the parameters of the model 
as well as by concentration of cars of different kind. This takes place if cars 
with higher maximum velocity are added to a system of slow cars. 

If all the cars have the smaller maximum velocity, I'max = 12, P2/P0 — 0.055, 
and n = 0.1, the probability of velocity has two maxima practically of the 
same height (Fig. 1). Changing these two parameters, the relative heights of 
the peaks change, and eventually one of them disappears, i.e., a first order 
phase transition is observed in the system. This phase transition is also driven 
by ratio of slow and fast cars (I'max = 20) as seen from Fig. 2 where the increase 
of concentration of fast cars to 30 % causes an increase of the small velocity 
peak with respect to the high velocity one. Further addition of the fast cars 
makes the system even slower and denser (Fig. 3). Only total disappearance 
of slow cars causes an abrupt increase of mean velocity and distances between 
the vehicles. Nevertheless, the flow is not free, because of large value of the 
parameter P2/P0 the system is in the congested phase (Fig. 4), with the most 
probable velocity much smaller than t'opt- For small values of the velocity 
spread P2/P0, potential /x, and large deceleration rate, the car flow freely, and 
the probability diagram consists of a peak at Vopt, narrow in F-direction and 
wide in X-direction. 

For the fast cars with maximum velocity ^^max = 20, the first order phase 
transition takes place at P2/P0 = 0.033, = 0.07 (Fig. 5). Now, addition of 
few slow cars, in distinction to the previous case, does not change the relative 
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Fig. 1. Probability density of headway X and car velocity V averaged over a group 
of 5 cars for a = 4, P2/V0 = 0.055, // = 0.1, and Umax = 12. 




16 

Fig. 2. Probability density of headway X and car velocity V averaged over groups 
of 5 cars for a = 4, P2/PQ = 0.055, = 0.1, with 70 % of cars of Umax = 12 and 30 % 

of Umax = 20. 
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Fig. 3. Probability density of headway X and car velocity V averaged over groups 
of 5 cars for a = 4, P2/P0 = 0.055, // = 0.1, with 2 % of cars of fmax = 12 and 98 % 

of ^^max = 20. 



Piy, 




Fig. 4. Probability density of headway X and car velocity V averaged over a group 
of 5 cars for o = 4, = 0.055, /x = 0.1, with 100 % of cars of Vmax = 20. 
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Fig. 5. Probability density of headway X and car velocity V averaged over a group 
of 5 cars for a = 4, P2/P0 = 0.033, // = 0.07, with 100 % of cars of t^max = 20. 

height of the two maxima, but totally destroys the probability diagram, and a 
wide peak in \^-direction and narrow in X-direction is formed. (The probabil- 
ity diagram is similar to that in Fig. 2.) Density is higher and velocity lower 
than of the system consisting only of slow cars (wmax = 12). 

As the cars cannot overtake each other, the most serious and instantaneous 
impact has a small admixture of cars with inferior properties to a free flowing 
pure system. On the other hand, addition of faster cars to a system of slow 
ones does not increase the mean velocity as well, but now the deterioration of 
the traffic flow is gradual. 

If the braking abilities of the cars are poor, they form platoons of vehicles with 
the same, low velocity. The density of cars in the platoons is not small, and 
is varying in a large interval. In Fig. 6, the deceleration rate a is 0.5 instead 
of a = 4 of previous cases. P2/P1 is 0.1, potential /i = 0.12, and the maximum 
velocity timax = 20. The presence of platoons is manifested by probability 
density maxima at integer values of the velocity. 

If few cars with good brakes (a = 4) are added, and the drivers are not able 
to distinguish the quality of the car ahead, they should be more careful and 
decrease their velocity to, approximately, one half. The drivers of the cars 
with poor brakes have to assume that all the cars have effective brakes and to 
increase the headways. The platoons disappear, and the increase of the mean 
velocity with the mean distance between cars is slow (Fig. 7). 
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Fig. 6. Probability density of headway X and car velocity V averaged over a group 
of 5 cars for Umax = 20, P2/P0 = 0.1, // = 0.12, with 100 % of cars of a = 0.5. 
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Fig. 7. Probability density of headway X and car velocity V averaged over groups 
5 indistinguishable cars for fmax = 20, P2/P0 = 0.1, /i = 0.12, with 90 % of cars of 
a = 0.5 and 10 % of a = 4. 
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Fig. 8. Probability density of headway X and car velocity V averaged over groups 
of 5 distinguishable cars for Vmax = 20, P2/P0 = 0.1, ^ = 0.12, with 90 % of cars of 
a = 0.5 and 10 % of a = 4. 

A different situation appears if the quality of the brakes of the preceding car 
is known to the driver. The cars with good brakes (type 1) may drive very 
closely to the cars with poor brakes (type 2). As the ratio of cars with good 
breaks is only 10% in Fig. 8, the cars with poor brakes drive mostly behind 
the cars of the same kind, the drivers know that, and formation of platoons 
is not hampered by the type 1 cars. All the cars move slowly so that also 
the distance between two cars with poor brakes is small, and the system is 
dense. Only the headway between type 1 and type 2 car is the same as all the 
headways in the previous case. 

For the same parameters, but a = 4, the cars move practically freely with 
only a small arm of slow vehicles (Fig. 9). This picture is destroyed by a small 
admixture (10%) of cars with poor brakes (Fig. 10). These cars should drive 
cautiously being surrounded by cars with good brakes which are able to stop 
immediately. They form very slow, dense platoons as at low velocity they are 
able to react effectively to slow cars in front of them. As in the front of cars 
with ineffective brakes there are practically always cars with good brakes, the 
possible knowledge of the type of the preceding car plays no role, and the 
diagram for P(V, X) remains the same for both above-mentioned cases. 

All the peaks corresponding to free and congested flow or formation of pla- 
toons result from correlations between the cars. It is seen, e.g., in Fig. 11, 
where the difference between the probability density D{y,X) of correlated 
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Fig. 9. Probability density of headway X and car velocity V averaged over a group 
of 5 cars for Umax = 20, P2/P0 = 0.1, /i = 0.12, with 100 % of cars of a = 4. 
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Fig. 10. Probability density of headway X and car velocity V averaged over groups 
of 5 cars for i^max = 20, P2/P0 = 0.1, /i = 0.12, with 10 % of cars of a = 0.5, 90 % 
of a = 4. 
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Fig. 11. Difference of probability densities (7) of correlated and uncorrelated group 
of 5 cars for the same parameters as in figure 5. 

and uncorrelated group of 5 cars, which P{V,X) is in Fig. 5, is shown. The 
probability density of a group of uncorrelated cars always displays only a single 
peak without any internal structure. 



4 Discussion 



Deriving canonical distribution of an equilibrium gas of classical particles in 
statistical mechanics, all the velocities and positions of the particles 
sumed to be equally probable, and a condition of given total mean energy is 
imposed on the system. (In grandcanonical distribution a condition for num- 
ber of particles is added.) In our approach to the gas of cars all their relative 
positions are equally probable as well, but conditional probabilities of their 
velocities are given by behaviour of drivers. Instead of the condition on total 
energy, equality condition of mean velocities of all cars and a condition on the 
mean total length of the system are imposed. Now, the probability of a state 
has not a simple form of an exponential of the Hamiltonian like in the case of 
physical systems. In principle, it can be rewritten into such a form, but now 
the quantities appearing in the exponential are not conserved in an isolated 
system, what is the case in a gas of classical particles in physics. That is why 
it is not possible to introduce an effective energy or Hamiltonian in our system 
of cars, i.e., our approach differs from that of thermodynamical traffic gas [12] 
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where interactions between cars and Hamiltonian of the system is introduced. 
Nevertheless, the results are similar but the velocity distribution of cars is not 
Gaussian. 

Our method is, to some extent, similar to a mean-field and cluster approxima- 
tion [8,9,10] looking for analytical solution of master equation corresponding 
to Nagel-Schreckenberg probabilistic cellular automaton model. Our choice 
of variables, velocities and headways, corresponds to that in the car oriented 
mean field theory [11], where the central role plays the probability Dn{v) iden- 
tical to our P{vi-i, Xi-i^i). In the above-mentioned approaches an approximate 
solution of master equations were found. In the present paper the solution is 
obtained by exact minimization of the entropy of the system. 

Our results are consistent with a number of other approaches, based on ki- 
netic or dynamical decripton of cars, and experimental observations of single 
line traffic [6,13,14]. It seems that systems of cars behave like many-particle 
physical systems or small parts of them, which, on microsopic level, are ruled 
by Newton or kinetic equations and, in equilibrium, are practically always 
described by canonical distribution derived from the principle of maximum 
entropy. 
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